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Abstract 

The field equations following from a Lagrangian L{R) will be de- 
duced and solved for special cases. If L is a non-linear function of 
the curvature scalar, then these equations are of fourth order in the 
metric. In the introduction we present the history of these equations 
beginning with the paper of H. Weyl from 1918, who first discussed 
them as alternative to Einstein's theory. In the third part, we give 
details about the cosmic no hair theorem, i.e., the details how within 
fourth order gravity with L = R + E? the inflationary phase of cosmic 
evolution turns out to be a transient attractor. Finally, the Bick- 
nell theorem, i.e. the conformal relation from fourth order gravity to 
scalar-tensor theory, will be shortly presented. 



*This paper is a siiortened version of my (still unpublished) "Lectures on Mathematical 
Cosmology" , see its preprint version gr-qc/0407095. Parts of my paper with R. Schimming 
"On the history of fourth order metric theories of gravitation" NTM-Schriftenr. Gesch. 
Naturw., Tech., Med. (Leipzig) 27 (1990) 41, gr-qc/0412038, are incorporated. 
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1 To the history of fourth order gravity 



Prom the advent of the general relativity theory (GRT) in 1915 by Albert 
Einstein (1879-1955) until today numerous geometrized theories of gravita- 
tion have been proposed. Here, we shall review the history of a class of 
theories which is conceptually rather close to GRT: 

- The gravitational field is described by a space-time metric only. 

- The field equation follows from a Hamiltonian principle. The Lagrangian 
L is a quadratic scalar in die Riemannian curvature of the metric. 

- The constants appearing in this ansatz are adjusted such that the theory 
is compatible with experimentally established facts. Hence, the Lagrange 
function reads 

L^aR^ + bRijR^ + kR + A (1.1) 

with constants a, b, k, A where a and b do not vanish simultaneously. The 
variational derivative of Rijki R^^^^ with respect to the metric can be linearly 
expressed by the variational derivatives of Rij R^^ and of R? [1] . Thus we 
may omit Riji^iR^^^^ in (1.1) without loss of generality. The theory is scale- 
invariant if and only if A ■ /c = 0. It is even conformally invariant if and only 
if A = A; = and 3a + 6 = 0. The field equation following from L eq. (1.1) 
is of fourth order, i.e. it contains derivatives up to the fourth order of the 
components of the metric with respect to the space-time coordinates. 

The fourth order metric theories of gravitation are a very natural modifi- 
cation of the GRT. Historically, they have been introduced as a specialization 
of Hermann Weyl's (1885-1955) nonintegrable relativity theory from 1918 [2]. 
Later on, just the fourth order theories became interesting and more and more 
physical motivations supported them: The fourth order terms can prevent 
the big bang singularity of GRT; the gravitational potential of a point mass 
is bounded in the linearized case; the infiationary cosmological model is a 
natural outcome of this theory. But all the arguments from classical physics 
were not so convincing as those from quantum physics: the quantization of 
matter fields with unquantized gravity background leads to a gravitational 
Lagrangian of the above form [3]. Moreover, fourth order theories turned 
out to be renormalizable at the one-loop quantum level [4], but at the price 
of losing the unitarity of the S-matrix. These circumstances caused a boom 



2 



of fourth order gravity (classical as well as quantum) in the seventies. We 
will stop our record of the history before this boom. We restrict ourselves 
to the purely metrical theories (i.e., the affinity is always presumed to be 
Levi-Civita) and want only to mention here that fourth order field equations 
following from a variational principle can be formulated in scalar-tensor the- 
ories, theories with independent affinity, and other theories alternative to 
GRT as well. 

1.1 Papers inspired by Weyl's theory 

In 1918, soon after Albert Einstein's proclamation of the GRT, Hermann 
Weyl proposed a new kind of geometry and a unified theory of gravitation 
and electromagnetism based on it. He dwelled on the matter in a series of 
papers [2, 5-10] until it became superseded by the modern gauge field inter- 
pretation of electromagnetism [11-13]. Note that the gauge concept together 
with the words "Eichung" (gauge) and "Eichinvarianz" (gauge invariancc) 
came into use in theoretical physics through Weyl's ansatz. For a broader 
discussion and evaluation we refer to [14]. A. Einstein [15] pointed out that 
the nonintegrability of the lengths of vectors under Weyl-like parallel prop- 
agation contradicts to physical experience. His argument has been refuted 
not earlier than in 1973: P. Dirac discusses the possibility of a varying grav- 
itational constant. He writes: 

"Such a variation would force one to modify Einstein's theory of grav- 
itation. It is proposed that the modification should consist in the revival 
of Weyl's geometry, in which lengths are nonintegrable when carried around 
closed loops, the lack of integrability being connected with the electromag- 
netic field" . [16, p. 403] 

H. Weyl's aesthetically very appeahng modification of GRT unfortunately 
does not directly describe the real dynamics of fields and particles; however it 
deeply infiuenced the "dynamics of theories" . By this we mean that various 
fundamental ideas have been formed or promoted by Weyl's papers: 

- the search for alternatives to the GRT based on geomctrization; 

- the unification of the interactions or forces of nature, beginning with gravity 
and electromagnetism; 
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- field theories based on tlie geometry of an affine connection; 

- conformal geometry and conformally invariant field theories; 

- the gauge field idea, and 

- fourth order gravitational field equations. 

Here we are interested just in the last item. Weyl required the La- 
grangian to be a polynomial function of the curvature and to be confor- 
mally invariant. He states: "Dies hat zur Folge, dass unscrc Thcoric wohl 
auf die Maxwellschen elektromagnetischen nicht abcr auf die Einstcinschen 
Gravitationsgleichungen fiihrt; an ihre Stelle treten Differentialgleichungen 
4. Ordnung." [2, p. 477] 

The ambiguity in the concrete choice of L appeared as a difficulty which 
is opposed to the spirit of unification: any linear combination of E? and 
Rij R^^ would do. The variation of Rij R^^ or of Rijki R^^^'' with respect to 
the vector field yields Maxwell-like equations, while for the choice of R^ an 
electromagnetic Lagrangian Fij F^^ together with a coupling constant a has 
to be added by hand: L — R^ + aFij F^^ [6, 2]. Weyl himself favoured differ- 
ent Lagrangians in different papers. Moreover, he took trouble to produce 
results compatible with Einstein's GRT. For this aim he destroyed the con- 
formal invariance by a special gauge. Ernst Reichenbacher criticizes: "Um 
so auffallender ist es, dass Weyl in dem von ihm durchgerechneten Beispiel 
fiir die Wirkungsfunktion durch Festlegung der Eichung vor der Variation 
den Grundsatz der Eichinvarianz durchbricht." [17, p. 157]. A more detailed 
analysis of the theory was necessary then. R. Weitzenbock [18] produced and 
studied all scalar invariants of the curvature in Weyl's geometry. W. Pauli 
[19, 20] and a little later F. Jiittner [21] calculated the spherically symmetric 
static gravitational field for variants of Weyl's theory. 

Pauli [20, S. 748] comes to an important conchision: "Hiernach ist klar, 
dass aus Beobachtungen der Merkurperihelbewcgung und der Strahlenablen- 
kung, die mit Einsteins Feldgleichungen im Einklang sind, niemals ein Ar- 
gument gegen Weyls Theorie entnommen werden kann, wenigstens solange 
die letztere eine der drei Invarianten R^ it!*-' , i?^, R^jki R^^^^ als Weltfunktion 
zugrunde legt." 

In other words, fourth order gravitational field equations following from 
(1.1) are not falsifiable by experimental physics! Pauli [20, 22] and other 
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authors did not even consider the vector field (i.e. assumed it to be equal to 
zero) in Weyl's theory, thus making it unafi^ected by the criticism of nonin- 
tegrabihty [15]. R. Bach [23] realized the possibihty to keep the conformal 
invariance in a purely metrical theory: a Lagrangian L — Cijki C"-''^' or equiv- 
alently, L — B? — 3Rij E}^ yields a conformally invariant field equation for 
the metric, later on called "Bach's equation" . In a similar spirit and in the 
same year 1921, Albert Einstein [24] proposed a conformally invariant theory. 
His expressions suffer from being non-rational in the metric. This theory is 
sometimes cited but has never been studied in details. 

Reichcnbachcr [25, 17] proposed a variant of Weyl's theory based on a 
non-rational Lagrangian resembling nonlinear Born-Infeld electrodynamics. 

In [26] , also L — is used to get a field-theoretical model for the electron, 
but the fourth order terms are lost by an error in the calculations. 

C. Lanczos [27] tried a programme of "Electromagnetismus als natiirhche 
Eigenschaft der Riemannschen Geometric" . He also assumed the vector field 
in Weyl's theory to be zero, but reintroduced it then in an alternative way as 
a set of Lagrangian multipliers. Unfortunately, Lanczos was, working with 
hyperbolic differential equations, misled by a formal analogy with elliptic 
differential equations. He varied the speculations with Lagrangian multipliers 
in a series of papers [28-34]. To take it positive, many useful mathematical 
formulas for fourth order theories resulted from Lanczos' work. Particularly, 
the paper [1] became a "citation classic" . 

In the twenties, the programme of classical field theory with its two cor- 
nerstones geometrization and unification lost some of its attractiveness in 
virtue of the quickly progressing quantum theory, cf. [35]. Moreover, there 
were the refutation of Weyl's theory and objections to fourth order equations. 
Lanczos expressed them as follows: 

"Der Grund, weshalb diese Untersuchungen nicht weiter gediehen sind 
und zu keinem Fortschritt fiihrten, lag an zwei Momenten. Einerseits war 
es entmutigend, dass man zumindest drei anscheinend gleichwertige Invari- 
anten zur Verfiigung hatte: B?, RapR"'^, R^p^sR"'^'^^ , ohne ein plausibles 
Auswahlprinzip zwischen ihnen zu besitzen. Andererseits erscheinen diese 
Gleichungen, solange man ihre innerc Struktur nicht verstehe, als Differen- 
tialgleichungen vierter Ordnung fiir die die Qik von einer Kompliziertheit sind. 



5 



die fiir jede weitere Schlussfolgerung ungeeignet ist." [27, p. 75] 

Similarly, Bergmann argues in his text-book [36] that, first, fourth order 
equations admit too many solutions and, second, their Lagrangian is rather 
ambiguous. This situation explains why only few papers on fourth order 
gravitation appeared in the period from the thirties to the sixties and why 
these did not follow the actual trends at their time. H. A. Buchdahl dealt with 
the subject in the period 1948-1980. In his papers he covered the following 
problems: 

- Invariant-theoretical considerations continuing those of Weitzenbock and 

Lanczos [37]; 

- General expressions for the variational derivatives of Lagrangians built from 
the curvature and, possibly, its derivatives are obtained [38-42]. 

- Einstein spaces are solutions of a rather general class of fourth order equa- 
tions [43, 44]; 

- Static gravitational fields in fourth order theories [45]; 

- Cosmological solutions in theories where the Lagrangian is a function of 

the scalar curvature [46, 47]; 

- Conformal gravity [48]; 

- Reinterpretation of some fourth order equations in five dimensions [49]. 

A. S. Eddington in 1921, see [50], and E. Schrodinger in 1948, see [51], also 
discussed gravitational field equations of fourth order to get field theoretical 
particle models, i.e., they tried to realize Einstein's particle programme. 

1.2 A new view 

Fourth order metric theories of gravitation have been discussed from 1918 
up to now. One original motivation was the scale invariance of the action, 
a property which does not hold in GRT. Another motivation was the search 
for a unification of gravity with electromagnetism, which is only partially 
achieved with the Einstein-Maxwell system. There was no experimental fact 
contradicting GRT which could give motives for replacing it by a more com- 
plicated theory. But a lot of problems appeared: 

1. The Lagrangian became ambiguous in sharp contrast to the required 
unification. 
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2. The higher order of the field equation brought mathematical problems 
(2.1.) in the search for solutions and physical problems (2.2.) for the inter- 
pretation of the additional degrees of freedom. 

3. The well-founded Newtonian theory of gravitation did not result as the 
weak-field limit of scale invariant fourth order gravity. 

The third problem was the last of these to be realized but the first to 
be solved, both in 1947: One has to break the scale invariancc of the theory 
by adding the Einstcin-Hilbert action to the purely quadratic Lagrangian. 
Then, up to an exponentially small term, the correct Newtonian limit appears 
[52]. 

The original scale invariant theory then, again emerges as the high-energy 
limit of that sum. The items 1., 2. and the absence of experimental facts 
contradicting GRT seemed to restrain die research on these theories already 
in the twenties. Only in 1966 a renewed interest in these theories arose 
in connection with a semiclassical description of quantum gravity [53-55]. 
The coefficient of the quadratic term became calculable by a renormalization 
procedure, thus solving problem 1, at least concerning the vacuum equa- 
tion. Further, the fact that fourth order gravity is one-loop renormalizable 
in contrast to GRT; a fact which was realized in 1977, [4] initiated a boom 
of research. It is interesting to observe that it is just the scale invariance of 
the curvature squared terms - the original motivation - which is the reason 
for the renormalizability. Also the latest fundamental theory - the super- 
string theory - gives in the field theoretical hmit (besides other terms) just 
a curvature-squared contribution to the action [56, 57]. The use of modern 
mathematics and computers has led to a lot of results to clarify the structure 
of the space of solutions thus solving problem 2.1. in the eighties. The more 
profound problem 2.2 has now three kinds of answers: 

a) In spite of the higher order of the differential equation, a prescribed 
matter distribution plus the 0(l/r)-behaviour of the gravitational poten- 
tial suffice - such as it takes place in Newtonian theory - to determine the 
gravitational potential for isolated bodies in a unique way for the weak-field 
slow-motion limit, [52, 53] 

b) the observation that the additional degrees of freedom are just the 
phases of damped oscillations which become undetectably small during the 
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cosmic evolution, and, by the way, can solve the missing mass problem and 
prevent the singularity problem of GRT [58] , and 

c) it is supposed that there exist massive gravitons besides the usual 
massless gravitons known from GRT, but they are very weakly coupled [59]. 

The last point to be mentioned is the experimental testability: In the 
recent three years many efforts have been made to increase the accuracy in 
determining the constants G, a and / if the gravitational potential is assumed 
(also in other theories than fourth order gravity) to be Gmr^^(\ + ae^''/') . 
The term proportional to a, the "fifth force" , can be interpreted as the fourth 
order correction to GRT. Up to now, it has not been possible to exclude a = 
by experiments [60-62]. 

In 1990, we concluded our paper on the history of fourth order gravity by 
saying: Fourth order gravity theories will remain an essential link between 
GRT and quantum gravity for a long time. Now, 16 years later, one can prove 
the correctness of this prediction by looking at references [63-99], [102], [105- 
109] which are all deal with a variant of fourth order gravity or its equivalents 
and its cosmological consequences and which is still far from being complete. 

2 Scalar fields and j{K) cosmology 

Following the paper "Comparing selfinteracting scalar fields and R + cos- 
mological models", gr-qc/0106035, see [82], we generalize the well-known 
analogies between m^0^ and R + theories to include the self-interaction 
A0'^-term for the scalar field. It turns out to be the R + R^ Lagrangian which 
gives an appropriate model for it. Considering a spatially fiat Friedmann 
cosmological model, common and different properties of these models are 
discussed, e.g., by linearizing around a ground state the masses of the corre- 
sponding spin 0-parts coincide. Then we prove a general conformal equiva- 
lence theorem between a Lagrangian L — L{R), 11 L" ^ 0, and a minimally 
coupled scalar field in a general potential. This theorem was independently 
deduced by several persons, and it is now known as Bicknell theorem [68] . 
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2.1 Introduction to scalar fields 



For the gravitational Lagrangian 

L = {R/2 + pR^)/87rG, (2.1) 
where /3 is some free but constant parameter, the value 

R=Rorit = -l/^P (2.2) 
is the critical value of the curvature scalar. It is defined by 

dL/dR = 0. (2.3) 



In regions where 
holds, we can define 
and 



R/Rcrit < 1 (2.4) 

Tj; = ln(l - R/Rcrit) (2.5) 



gij = {l- R/ i?crit) 9ij ■ (2.6) 

In units where SttG = 1 we now obtain from the Lagrangian eq. (2.1) via 
the conformal transformation eq. (2.6) the transformed Lagrangian 

L^R/2- 'ifH-^.il^ - (l - e-^)' /16/3 (2.7) 

being equivalent to L, cf. Whitt [85]. 

For /3 < 0, i.e., the absence of tachyons in L eq. (2.1), we have massive 
gravitons of mass 

mo = (-12/?)-i/^ (2.8) 

in L, cf. Stelle [4]. For the weak field limit, the potential in eq. (2.7) can 
be simplified to be t/'^/(16 • /3), i.e., we have got a minimally coupled scalar 
field whose mass is also mo. The superfiuous factor 3/2 in eq. (2.7) can 
be absorbed by a redefinition of ■0- Therefore, it is not astonishing, that 
all results concerning the weak field limit for both R + i?^-gravity without 
tachyons and Einstein gravity with a minimally coupled massive scalar field 
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exactly coincide. Of course, one cannot expect this coincidence to hold for 
the non-linear region, too, but it is interesting to observe which properties 
hold there also. 

We give only one example here: we consider a cosmological model of the 
spatially flat Priedmann type, start integrating at the quantum boundary, 
which is obtained by Rij^iR!'^^^ on the one hand, and Tqo on the other hand, 
to have Planckian values, with uniformly distributed initial conditions and 
look whether or not an inflationary phase of the expansion appears. In both 
cases we get the following result: The probability p to have sufficient inflation 
is about p = 1 — V^^^o/''^Pb i-G-; P = 99.992% if wc take rriQ = lO^^mpi from 
GUT and A = 64, where is the linear multiplication factor of inflation. 

Prom Quantum field theory, however, instead of the massive scalar field, 
a Higgs field with self-interaction turns out to he a better candidate for 
describing effects of the early universe. One of the advances of the latter 
is its possibility to describe a spontaneous breakdown of symmetry. In the 
following, we try to look for a purely geometric model for this Higgs field 
which is analogous to the above mentioned type where L — R + modelled 
a massive scalar fleld. 



The ground states are deflned by ^ = const, and dL/d(j) = 0. This means 
= for the scalar field, and the three ground states = 0o = 0, and 



2.2 The Higgs field 



Por the massive scalar Pield (f) we have the mater Lagrangian 




(2.9) 



and for the Higgs field to be discussed now. 




(2.10) 




(2.11) 



for the Higgs field. 
The expression 




(2.12) 
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represents the effective mass at these points. This gives the value m for the 
scalar field eq. (2.9), so justifying the notation. Further, eq. (2.12) give mass 
i /i at = and \/2/i at = 0± for the Higgs field eq. (2.10). The imaginary 
value of the mass at the ground state = shows the instability met there, 
and in the particle picture, this gives rise to a tachyon. 

To have a vanishing Lagrangian at the ground state 0± eq. (2.11) we add 
a constant 

A = -3//V2A (2.13) 
to the Lagrangian eq. (2.10). The final Lagrangian reads 

L = R/2 + Lx + A (2.14) 

with La eq. (2.10) and A eq. (2.13). 

2.3 The non-linear gravitational Lagrangian 

Preliminarily we direct the attention to the following fact: on the one hand, 
for Lagrangians (2.9), (2.10) and (2.14) the transformation (p — > —0 is a pure 
gauge transformation, it does not change any invariant or geometric objects. 
On the other hand, 

Rijkl —Rijki (2-15) 

or simpler 

R^-R (2.16) 
is a gauge transformation at the linearized level only: taking 

Qik = Vik + ehik , (2.17) 

where 

7/,fc = diag(l, -1, -1, -1) , (2.18) 

then e — >■ —e implies curvature inversion eq. (2.15). To be strict: this is valid 
at the linearized level in e only. On the other hand, the curvature inversion 
eq. (2.15), and even its simpler version eq. (2.16), fails to hold quadratic in 
e. This corresponds to fact that the £^-term in eq. (2.7), which corresponds 
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to the ^^-term in the development of L in powers oi ip, is the first one to 
break the ■?/' — > —ip symmetry in eq. (2.7). 
In fact, the potential is essentially 

(l-e(-^^)\ (2.19) 
Calculating this, the result reads 

+ + (2.20) 

Now, let us introduce the general non-linear Lagrangian L = L{R) which 
we at the moment only assume to be an analytical function of R. The ground 
states are defined by i? = const., i.e., 

L'Rik-9ikL/2^0. (2.21) 

Here, L' = dL/dR. 



2.4 Calculation of the ground states 

From eq. (2.21) one immediately sees that dL/dR = defines critical values 
of the curvature scalar. For these values R = Rait it holds: For L(i?crit) 7^ 
no such ground state exists, and for L(i?crit) = 0, we have only one equation 
R — -Rcrit to be solved with 10 arbitrary functions Qik- We call these ground 
states degenerated ones. For L — R^, Rcru — 0, this has been discussed by 
Buchdahl [41]. Now, let us concentrate on the case dL/dR ^ 0. Then i?^ is 
proportional to gij with a constant proportionality factor, i.e., each ground 
state is an Einstein space 

R,, = Rg,,/^, (2.22) 

with a prescribed constant value R. Inserting eq. (2.22) into eq. (2.21) we 
get as condition for ground states 

RU = 2L . (2.23) 

As an example, let L be a third order polynomial 

L = A + i?/2 + + Ai?Vl2 . (2.24) 
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Wc consider only Lagrangians with a positive linear term as we wish to 
reestablish Einstein gravity in the A — > weak field limit, and f3 < to 
exclude tachyons there. 

We now solve eq. (2.23) for the Lagrangian eq. (2.24). For A = we have, 
independently of /3, the only ground state R — — 4A. It is a degenerated one 
if and only if /3A = 1/16. That implies that for usual R -\- B? gravity eq. 
(2.1), i.e. A = A = 0, we get R — Q as the only ground state; it is a 
non-degenerated one. 

Now, let A 7^ and A = 0. To get non-trivial ground states wc need the 
additional assumption A > 0. Then, besides i? = 0, the ground states are 



being quite analogous to the ground states eq. (2.11) of the Higgs field 
eq. (2.10). The ground state i? = is not degenerated. Of course, this 
statement is independent of A and holds true, as one knows, for A = 0. To 
exclude tachyons, we require /? < 0, then i?_ is not degenerated and i?+ 
is degenerated if and only li (5 — The case AA will not be 

considered here. 

2.5 The field equation 

For L = L[R) the variation 



gives with L' = OL/dR the following fourth order gravitational field equation 




(2.25) 




(2.26) 



L'Ri^ - + g,pL' - L'.,^ = . 



(2.27) 



It holds 



L'jj — L"R-ij + L'" R-iR-j ■ 
With eq. (2.28), the trace of eq. (2.27) reads 



(2.28) 



L'R -2L + 3L"nR + 3L'" R.^K'^ = , 



(2.29) 
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i.e., with L eq. (2.24) 

-2A - R/2 + \R^/12 + &(3UR + —{RUR + R-kR^) = . (2.30) 

Inserting eqs. (2.32), (2.33) and (2.28) into the 00-component of eq. (2.27) 
we get the equation 

= h^l2 - A/6 - 6/3(2/i/i - /i^ + 6/i^/i) 
+3A(/i + 2h^) {6hh + 19h^h - 2h^ - 2/i^) . (2.31) 

The remaining components are a consequence of this one. 

2.6 The Friedmann model 

Now we take as Lagrangian eq. (2.24) and as hne element 

(is2 = (1^2 _ a2(^^^(^^^2 ^ ^ ^^2^) _ (^2.32) 

The dot denotes d/dt and h — d/a and we have 

R = -Qh - 12h^ . (2.33) 

Here we only consider the spatially flat Friedmann model eq. (2.32). There- 
fore, we can discuss only de Sitter stages with R < 0, especially the ground 
state i?+ eq. (2.25) representing an anti-de Sitter spacetime does not enter 
our discussion, but -R_ does. 

Now, let A = 0. Solutions of eq. (2.31) with constant values h are h = 
representing fiat spacetime and in the case that A > also 

h = (2.34) 

representing the dc Sitter spacetime. These are the non-degenerated ground 
states = and R = R- = — \/6/A, respectively. Eq. (2.31) can be written 
as 

= h\l - 24A/i^) /2 + hh (l/ml + 18X{h + 2h^)) 
-6Xh^ + h'^i^bXh? - I /2ml) + 3/i^/i(l/m^ + ZQXh'^) . (2.35) 
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First, let us consider the singular curve defined by the vanishing of the coef- 
ficient of h in eq. (2.35) in the h — /i-phase plane. It is, besides /i = 0, the 
curve 

-2h^ -l/l%\ml (2.36) 

i.e., just the curve 

i?= l/3Am^ = -4/?/A (2.37) 

which is defined by L" = 0, cf. eq. (2.29). This value coincides with i?+ if 
P — — y^3A/8, this value we need not discuss here. Points of the curve eq. 
(2.36) fulfil eq. (2.35) for 

h = ±l/ ISXmlVs \Jl-l jlSXml (2.38) 

only, which is not real because of A <^ rriQ. Therefore, the space of solutions 
is composed of at least two connected components. 
Second, for /i = we have /i = or 

/i = -l/12Am^ (2.39) 

Prom the field equation we get: h — h — implies hh > 0, i.e. h does not 
change its sign. In a neighbourhood of eq. (2.39) we can make the ansatz 

oo 

h = -t/12Xml + ^ a„ r (2.40) 

n=2 

which has solutions with arbitrary values 02. This means: one can change 
from expansion to subsequent recontr action, but only through the "eye of 
a needle" eq. (2.39). On the other hand, a local minimum of the scale 
factor never appears. Purther, eq. (2.39) does not belong to the connected 
component of flat spacetime. 

But we are especially interested in the latter one, and therefore, we restrict 
to the subset h > h{ cq. (2.36)) and need only to discuss expanding solutions 
h > 0. Inserting h = 0, 

h = /i(24A/i^ - I) /{2/ml + 72A/i2) (2.41) 

turns out, i.e., > for /i > l/v^24A only. All other points in the h — h phase 
plane are regular ones, and one can write dh/dh = h/h — F{h,h) which can 
be calculated by eq. (2.35). 
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For a concrete discussion let A ~ lO^/pj and mo = 10^'''mpi. Then both 
conditions j3 <^ —\/~\ and \R-\ < Ipi are fulfilled. Now the qualitative be- 
haviour of the solutions can be summarized: There exist two special solutions 
which approximate the ground state R- for t — > — oo. All other solutions 
have a past singularity h ^ oo. Two other special solutions approximate the 
ground state i?_ for t — > +00. Further solutions have a future singularity 
h 00, and all other solutions have a power-like behaviour for t 00, a{t) 
oscillates around the classical dust model a{t) ~ t^^^. But if we restrict the 
initial conditions to lie in a small neighbourhood of the unstable ground state 
-R_, only one of the following three cases appears: 

1. Immediately one goes with increasing values /i to a singularity. 

2. As a special case: one goes back to the de Sitter stage 

3. The only interesting one: One starts with a finite Zpi-valued infiationary 
era, goes over to a GUT-valued second infiation and ends with a power-hke 
Friedmann behaviour. 

In the last case to be considered here, let A = 0, A > and (3 < 0. The 
analogue to eq. (2.35) then reads 

= /1V2 - A/6 + {2hh -h^ + eh^ii) /2ml ■ (2-42) 

Here, always h ^ holds, we consider only expanding solutions h > 0. For 
h — we have 

h = {Aml/3 - mlh^) /2h . (2.43) 
For h — we have h > ml/G and 

h = (A/3 + hymiy/\l + 6/i/m2)-i/2 . (2.44) 

We obtain the following result: All solutions approach the de Sitter phase 
= A/3 as t — 00. There exists one special solution approaching h = 
—771^/6 for h 00, and all solutions have a past singularity h 00. For 
a sufficiently small value A we have again two different inflationary eras in 
most of all models. 
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2.7 The generalized Bicknell theorem 



In this section wc derive a general equivalence theorem between a non-linear 
Lagrangian L{R) and a minimally coupled scalar field with a general po- 
tential with Einstein's theory. Instead of we take 

^ = -^2/3 0. (2.45) 

This is done to avoid square roots in the exponents. Then the Lagrangian 
for the scalar field reads 



L^R/2- 3g'^iP,iip.j/A + ViiP) . (2.46) 
At ground states ip — ipo, defined by dV/dip = the effective mass is 

m ^ y/2/3^/¥v/di^ , (2.47) 
cf. eqs. (2.12) and (2.45). The variation = SL/5ip gives 

= dV/dij + 3 (2.48) 
and Einstein's equation is 

Eij = KTij (2.49) 

with 

KTij = 3V';^V';,/2 + Qij (v{iIj) - -r''M;b) . (2.50) 

Now, let 



~9ij^e^9ij- (2.51) 

The conformal transformation eq. (2.51) shall be inserted into eqs. (2.48), 
(2.49) and eq. (2.50). One obtains from eq. (2.48) with 

:= g'''^.i 

n^jj + ^'%k = -2{e'f'dV/d^jj) /3 (2.52) 
and from eqs. (2.49), (2.50) 

Eij = + + Qij (e'^1/(V') - Ot/j - V^.^^'") ■ (2.53) 
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Its trace reads 

-R = Ae^Viij) - SDijj - 3^;„^'" . (2.54) 
Comparing with eq. (2.52) one obtains 

R = R{^) = -2e-^9 (e'^F(^)) /Ot/j . (2.55) 

Now, let us presume dR/dip ^ 0, then eq. (2.55) can be inverted as 

= F{R) . (2.56) 

In the last step, eq. (2.56) shall be inserted into eqs. (2.52), (2.53), (2.54). 
Because of 

F{R).ij = dF/dR ■ R,ij + d'^F/dR^ ■ R,iR,j (2.57) 

and dF/dR ^ 0, eq. (2.53) is a fourth-order equation for the metric Qij. We 
try to find a Lagrangian L = L{R) such that the equation 5L^J—g/5g^^ = 
becomes just eq. (2.53). For L' = OL/dR ^ 0, eq. (2.27) can be solved to 
be 

Eij = -9^jR/2 + g^jL/2L' - g,pL' / L' - L'.,^/L' . (2.58) 

We compare the coefficients of the R.^ij terms in eqs. (2.53) and (2.58), this 
gives 

dF/dR = L"/L' , hence 

LiR) r e^^^'^dx + Ao (2.59) 
JRo 

with suitable constants Aq, and Rq, /j, ^ 0. We fix them as follows: We 
are interested in a neighbourhood of R = Rq and require L'[Rq) — 1/2. 
Otherwise L should be multiplied by a constant factor. Further, a constant 
translation of if) can be used to obtain F{Rq) = 0, hence = 1/2, L{Rq) = 
Aq, and 

L'{Rq) = dF/dR{RQ)/2 ^ . (2.60) 

With eq. (2.59) being fulfilled, the traceless parts of eqs. (2.53) and (2.59) 
identically coincide. Furthermore, we have 

UL'/L' = DF + F''F.i (2.61) 
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and it suffices to test the validity of the relation 

V(F(R)) = -R/2 + L/2L' . (2.62) 

It holds 

2L'^e^, i.e., 

e'-^l/(F(i?)) -L-i?e-^/2. (2.63) 

At R = Ro, this relation reads V{0) = Aq — i?o/2- Applying d/dR to eq. 
(2.63) gives just eq. (2.52), and, by the way, V'{0) — Rq/2 — 2Kq. In sum, 

LiR) = VU]) + Ro/2 + r e^(^)rfa;/2 , (2.64) 

JRo 

where F{x) is defined via F{Rq) = 0, 

= F (-2e-'^d{e''''V{^l^))/d^f) . (2.65) 

Now, let us go the other direction: Let L — L{R) be given such that at 
R = i?o, L'L" 7^ 0. By a constant change of L let L'(i?o) = 1/2. Define 
Aq = L{Rq), — F{R) — ln(2L'(i?)) and consider the inverted function 
i? = F-i(V'). Then 

Vm = (Ao - Ro/2)e-^'^ - e"'^ T F-Hx)dx/2 (2.66) 

is the potential ensuring the above mentioned conformal equivalence. This 
procedure is possible at all i?-intervals where V L" ^ holds. For analytical 
functions L{R), this inequality can be violated for discrete values R only, or 
one has simply a linear function L{R) being Einstein gravity with A-term. 

It turned out that this integral eq. (2.66) can be evaluated in closed form 
as follows: 

Vii/j) = L{F-\ij)) ■ e-'^ - - F-\^) ■ e-^ . (2.67) 
Examples: 1. Let L = A + , Rq = 1/4, then AR = e'^ and 

y(^) = Ae-2'^-l/16. (2.68) 
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2. Let L = A + RI2 + I3R^ + \R^/12, Rq = 0, hence /3 7^ is necessary. 
We get 

- 1 = A(3R + \R^/2 and 

2/3A-^e-2^ (e^ - 1 - 16/3'(3A)-^((1 + A(e^ - I) /8/3')=^/2 _ _ (2.69) 

The hmit A ^ in eq. (2.69) is possible and leads to 

V{iP) = Ae-'^ - (1 - e-'^) VI6/3 , (2.70) 

so we get for A = again the potential from eq. (2.7). 
Now, let Rq be a non-degenerated ground state, hence 

L{R) = Ao + (i? - i?o)/2 + L"(i?o)(i? - i?o)V2 + ■ ■ ■ (2.71) 

with L"{Rq) 7^ and Aq = -Ro/4, cf. subsection 2.4. Using eq. (2.66) we get 
V'{Q) = and 

F"(0) = i?o/2-l/(4L"(i?o))- (2.72) 
2.8 The generalized equivalence 

It is worth mentioning that this conformal equivalence theorem can be for- 
mulated for arbitrary dimensions n> 2: Let 

C = R/2- \g'^\^4>\j + V{<P) (2.73) 

and 

Qii = e^^ga , A = (2.74) 

(n-l)(n-2) ^ ^ 

be the conformally transformed metric. Then the solutions of the variation 
of eq. (2.73) are transformed by eq. (2.74) to the solutions of the variation 
of L = L{R) , where 
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is supposed to be locally, i.e. near R = Rq, invertible as 

4>^F{R), FiR^)^0, F'{Ro)y^O 

C{R) = Iro + V{0) + 1 r e^^'^^/'^dx . (2.76) 
2 2 Jro 

The inverse direction is possible provided L'{R) ■ L"{R) 7^ 0. 
2.9 The fourth order gravity model 

Now we come to the following question: Let L[R) be a C^-function fulfilling 
L(0) = 0, L'(0) • L"(0) < 0. Then we can write 

C{R) ^^ + I3R^ + 0{R^) , P<0. (2.77) 

We consider the Bianchi-type I vacuum solutions which start in a neigh- 
bourhood of the Minkowski spacetime and ask for the behaviour as t — > 00. 
Applying the equivalence theorem we arrive at the models discussed before, 
and this is applicable for \R \ being small enough. The conformal factor de- 
pends on t only, and therefore, the space of Bianchi-type I models will not be 
leaved, and we can formulate the following: In a neighbourhood of Minkowski 
spacetime, all Bianchi-type I models which represent a stationary point of 
the action eq. (2.77), can be integrated up to t ^ 00 or —00, let it be +00. 
One singular solution is the Kasncr solution and all other solutions undergo 
isotropization and have an averaged equation of state p = for t — > 00. 
The de Sitter spacetime 

ds' = dt^ - e'^*((ia;' + dy^ + dz^) , H^O (2.78) 

is the spacetime being mainly discussed to represent the inflationary phase 
of cosmic evolution. However, a spacetime defined by 

ds^ ^dt^ -\t\^P{dx^ + dy^ + dz^), py^O (2.79) 

enjoys increasing interest for these discussions, too. EspeciaUy, eq. (2.79) 
with p > 1, i > is called power-law infiation; and with p < 0, i < it is 
called polar infiation. 
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Eq. (2.79) defines a self-similar spacetime: if we multiply the metric 
ds^ by an arbitrary positive constant a^, then the resulting ds"^ = o?ds^ is 
isometric to ds^. Proof: We perform a coordinate transformation i — at, 
X — b{a,p)x . . . On the other hand, the de Sitter spacetime eq. (2.78) is not 
self-similar, because it has a constant non-vanishing curvature scalar. 

Power-law inflation is intrinsically time-oriented. Proof: The gradient of 
the curvature scalar defines a temporal orientation. On the other hand, the 
expanding {H > 0) and the contracting {H < 0) dc Sitter spacetime can be 
transformed into each other by a coordinate transformation, because both of 
them can be transformed to the closed Friedmann universe with scale factor 
cosh(ift), which is an even function of t. This property is connected with 
the fact that eq. (2.79) gives a global description, whereas eq. (2.78) gives 
only a proper subset of the full de Sitter spacetime. 

For p — > oo, eq. (2.79) tends to eq. (2.78). Such a statement has to be 
taken with care, even for the case with real functions. Even more carcftilly 
one has to deal with spacetimes. The most often used limit - the Gcroch- 
limit of spacetimes - has the property that a symmetry, here we take it to be 
self-similarity, of all the elements of the sequence must also be a symmetry 
of the limit. 

From this it follows that the Geroch limit of spacetimes eq. (2.79) with 
p — > oo cannot be unique, moreover, it is just the one-parameter set eq. 
(2.78) parametrized by arbitrary values H > 0. 

3 Cosmic no hair theorem and Newtonian 
limit 

We now discuss the cosmic no hair theorem, which tells under which circum- 
stances the de Sitter spacetime represents an attractor solution within the set 
of other nearby solutions. This property ensures the inflationary model to be 
a typical solution. The notion "cosmic no hair theorem" is chosen because 
of its analogous properties to the "no hair theorem" for black holes. 

After a general introduction we restrict our consideration to spatially flat 
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Friedmann models. In this section, we choose gravitational Lagrangians 



and linear combinations of them. They are motivated from trials how to 
overcome the non-renormalizability [3, 4] of Einstein's theory of gravity. 
Results are: For arbitrary k, i.e., for arbitrarily large order 2A; + 4 of the field 
equation, one can always find examples where the attractor property takes 
place. Such examples necessarily need a non- vanishing i?^-term. The main 
formulas do not depend on the dimension, so one gets similar results also for 
1+1-dimensional gravity and for Kaluza-Klein cosmology. 

Over the years, the notion "no hair conjecture" drifted to "no hair theo- 
rem" without possessing a generally accepted formulation or even a complete 
proof. Several trials have been made to formulate and prove it at least for 
certain special cases. They all have the overall structure: "For a geometri- 
cally defined class of spacetimes and physically motivated properties of the 
energy-momentum tensor, all the solutions of the gravitational field equation 
asymptotically converge to a space of constant curvature." 

Weyl wrote: "The behaviour of every world satisfying certain natural 
homogeneity conditions in the large follows the de Sitter solution asymp- 
totically." Barrow and Gotz apply the formulation: "All ever-expanding 
universes with A > approach the dc Sitter spacetime locally." 

Barrow gave examples that the no hair conjecture fails if the energy con- 
dition is relaxed and pointed out, that this is necessary to solve the graceful 
exit problem. He uses the formulation of the no hair conjecture "in the pres- 
ence of an effective cosmological constant, stemming e.g. from viscosity, the 
de Sitter spacetime is a stable asymptotic solution" . 

In the three papers [80], Prigoginc ct al. developed a phenomenological 
model of particle and entropy creation. It allows particle creation from space- 
time curvature, but the inverse procedure, i.e. particle decay into spacetime 
curvature is forbidden. This breaks the t — > — t-invariance of the model. 
Within that model, the expanding de Sitter spacetime is an attractor solu- 
tion independently of the initial fiuctuations; this means, only the expanding 
de Sitter solution is thermodynamically possible. 

Vilenkin discussed future-eternal infiating universe models: they must 




(3.1) 
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have a singularity if the condition D: "There is at least one point p such that 
for some point q to the future of p the volume of the difference of the pasts 
of p and q is finite" is fulfilled. 

In Maeda [77] the following argument is given: If the matter distribution 
is too clumpy, then a large number of small black holes appears. Then one 
should look for an inflationary scenario where these black holes are harmless. 
They cannot clump together to one giant black hole because of the exponen- 
tial expansion of the universe; this explains the existing upper bound of the 
mass of black holes in the quasi- de Sitter model: above 

Merit = (3.2) 

there do not exist horizons; this restriction is called cosmic hoop conjecture. 

For the no hair conjecture for E? models one uses the formulation "asymp- 
totical de Sitter, at least on patch" . The restriction "on patch" is not strictly 
defined but refers to a kind of local validity of the statement, e.g., in a region 
being covered by one single synchronized system of reference in which the 
spatial curvature is non-positive and the energy conditions are fulfilled. The 
Starobinsky model [103] is one of those which does not need an additional 
inflaton field to get the desired quasi de Sitter stage, see [100, 101, 104] as 
three of the earliest approaches to this point. 

A further result is that by the addition of a cosmological term, the 
Starobinsky model leads naturally to double inflation. Let us comment this 
result: It is correct, but one should add that this is got at the price of getting 
a "graceful exit problem" , by this phrase there is meant the problem of how 
to finish the inflationary phase dynamically - in the Starobinsky model this 
problem is automatically solved by the fact that the quasi de Sitter phase is 
a transient attractor only. 

The paper Buchdahl [38] deals with Lagrangians of arbitrarily high order. 
Its results are applied in [63, 65, 81] to general Lagrangians □). The 

original idea of A. D. Sakharov [54] from 1967 was to define higher order 
curvature corrections to the Einstein action to get a kind of elasticity of the 
vacuum. Then the usual breakdown of measurements at the Planck length, 
such a short de Broglie wave length corresponds to such a large mass which 
makes the measuring apparatus to a black hole, is softened. 
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3.1 Definition of an asymptotic de Sitter spacetime 

In this subsection we want to compare some possible definitions of an asymp- 
totic de Sitter spacetime. To tliis end let us consider the metric 

n 

ds^^df-e^'^^'^J^di^'f (3-3) 
1=1 

which is the spatially fiat Friedmann model in n spatial dimensions. We 
consider all values n > 1, but then concentrate on the usual case n = 3. The 
Hubble parameter is 

H^a = ^. (3.4) 
dt ^ ^ 

We get 

Roo = -n + H^^ , R = -2nH - n{n + \)E^ . (3.5) 

Then it holds: The following 4 conditions for metric (3.3) are equivalent. 
1: The spacetime is flat. 2: It holds R — Rqo — 0. 3: The curvature 
invariant RijR^^ vanishes. 4: Either a — const, or 

n — 1 and a = In |i — io| + const. (3.6) 

For the last case with n—1 one has to observe that 

ds^ = dt^ -{t- tofdx^ (3.7) 

represents flat spacetime in polar coordinates. For the proof we use the 
identity 

R,jR^ = {R^o? + -{R- RooY ■ (3.8) 
n 

An analogous characterization is valid for the de Sitter spacetime. The 
following 4 conditions for metric (3.3) arc equivalent. 1: It is a non-flat 
spacetime of constant curvature. 2: i?oo = R/ {n + 1) = const. 7^ 0. 3: 
(n + l)RijR'^ = R^ = const, 0. 4: Either H = const, or 

n = 1 and (ds^ = dt^ - sm^{Xt)dx^ or ds^ = dt^ - sinh^ {Xt)dx^) . (3.9) 
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For n = 1, the de Sitter spacetime and anti-de Sitter spacetime differ by the 
factor — 1 in front of the metric only. For n > 1, only the de Sitter spacetime, 
having i? < 0, is covered, because the anti-de Sitter spacetime cannot be 
represented as spatially fiat Friedmann model. Our result shows that within 
the class of spatially flat Friedmann models, a characterization of the de 
Sitter spacetime using polynomial curvature invariants only, is possible. 

Next, let us look for isometrics leaving the form of the metric (3.3) in- 
variant. Besides spatial isometrics, the map a — > a defined by 

a{t) = c + a{±t + to), c, to — const. (3.10) 

leads to an isometric spacetime. The simplest expressions being invariant by 
such a transformation are and H. We take a as dimensionless, then H is 
an inverse time and H an inverse time squared. Let H ^ in the following. 
The expression 

£ := HH-^ (3.11) 

is the simplest dimensionless quantity defined for the spatially fiat Friedmann 
models and being invariant with respect to this transformation. Let n > 1 
in the following, then it holds: Two metrics of type (3.3) are isometric if and 
only if the corresponding functions a and a are related by this transforma- 
tion. It follows: Metric (3.3) with H represents the de Sitter spacetime 
if and only if £ = 0. 

All dimensionless invariants containing at most second order derivatives 
of the metric can be expressed as f{s), where / is any given function. But if 
one has no restriction to the order, one gets a sequence of further invariants 

e2 = HH-\ e^=-^H-^-K (3.12) 

Let if > in metric (3.3) with n > 1. We call it an asymptotic de Sitter 
spacetime if 

lim ^ = const. ^ (3.13) 



or 



lim = const. > and hm {n + l)RijRij - R^ ^ (3.14) 

t—^oo t— »oo 
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or for some natural number p one has 

lim = for l<j<p. (3.15) 

In general, all these definitions are different. Using the identity 

RijR'^ = n^(H + H'^f + n{H + nH'^f (3.16) 

we will see that all these definitions lead to the same result if we restrict our- 
selves to the set of solutions of the higher-order field equations. The problem 
is that none of the above definitions can be generalized to inhomogeneous 
models. One should find a polynomial curvature invariant which equals a 
positive constant if and only if the spacetime is locally the de Sitter space- 
time. To our knowledge, such an invariant cannot be found in the literature, 
but also the non-existence of such an invariant has not been proven up to 
now. 

This situation is quite different for the positive definite case: For signature 
(+ + ++) and Sij = Rij - jgij it holds: 

I={R- Rof + QjkiC'^'' + SijS'^ = (3.17) 

if and only if the V4 is a space of constant curvature Rq. So I — > is a 
suitable definition of an asymptotic space of constant curvature. One possi- 
bility exists, however, for the Lorentz signature case, if one allows additional 
structure as follows: An ideal fluid has an energy-momentum tensor 

Tij = {p + p)uiUj - pQij (3.18) 

where Ui is a continuous vector field with UiU^ = 1. For matter with equation 
of state p = —p. the equation T^-'.j = implies p = const., and so every 
solution of Einstein's theory with such matter is isometric to a vacuum so- 
lution of Einstein's theory with a cosmological term. The inverse statement, 
however, is valid only locally: 

Given a vacuum solution of Einstein's theory with a A-term, one has 
to find continuous timelike unit vector fields which need not to exist from 
topological reasons. And if they exist, they are not at all unique. So, it 
becomes possible to define an invariant J which vanishes if and only if the 
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spacetime is de Sitter by transvccting the curvature tensor with u^u^ and/or 
g''^ and suitable hnear and quadratic combinations of such terms. Then time 
t becomes defined by the streamhnes of the vector u^. If one defines the 
asymptotic de Sitter spacetime by J — > as i — > oo, then it turns out, 
that this definition is not independent of the vector field u\ 

3.2 Lagrangian F(i?, Di?, D^i?, . . . , n^R) 

Let us consider the Lagrangian density L given by 

L = F{R, DR, n^R, n''R)yf^ (3.19) 

where R is the curvature scalar, □ the D'Alembertian and gij the metric 
of a Pseudo-Riemannian Vr, of dimension D > 2 and arbitrary signature; 
g — — \ det gij \ . The main application will be D = 4 and metric signature 
(H ). F is supposed to be a sufficiently smooth function of its argu- 
ments, preferably a polynomial. Buchdahl already dealt with such kind of 
Lagrangians in 1951, but then it became quiet of them for decades. 

The variational derivative of L with respect to the metric yields the tensor 

P'^ = ^ ^ (3.20) 

The components of this tensor read 

Pij — GRij — -jFgij — G-ij + gijOG + Xij (3-21) 

where the semicolon denotes the covariant derivative, Rij the Ricci tensor, 
and 

X^J = E l9iAFA{°^-'R)n,m - FA^.,[n^-'R].,) (3.22) 

A=l ^ 

having the round symmetrization brackets in its last term. For /c = 0, i.e. 
F — F{R), a case considered in subsection 3.3, the tensor X^j identically 
vanishes. It remains to define the expressions Fa, A — 0, . . . ,k . 
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and for ^4 = A; — 1, . . . , 

dF 

Fa = OFa+1 + (3.24) 

and finally G — Fq. The brackets are essential, for any scalar $ it holds 

□ - = i?/' (3.25) 

Inserting $ = D™/? into this equation, one gets identities to be applied in 
the sequel without further notice. The covariant form of energy-momentum 
conservation reads 

P% = (3.26) 

and Pij identically vanishes if and only if F is a divergence, i.e., locally there 
can be found a vector such that F — v\, holds. Remark: Even for com- 
pact manifolds without boundary the restriction "locally" is unavoidable, for 
example, let D = 2 and V2 be the Riemannian two-sphere S'^ with arbitrary 
positive definite metric. Ris a divergence, but there do not exist continuous 
vector fields fulfilling i? = on the whole 5"^. Example: form, n > 
it holds 

0"'R □"i? - R = divergence . (3.27) 

So, the terms 0"^R with naturals m and n can be restricted to the case 
m — without loss of generality. 



3.3 No hair theorems for higher-order gravity 

For n > 1, the n + 1-dimensional de Sitter spacetime is an attractor solution 
for the field equation derived from the Lagrangian 

R(^+^y^ . (3.28) 

It is not an attractor solution for the Lagrangian Rai^R and A: > 0. There 
exist combinations of coefficients q, such that the de Sitter spacetime is 
an attractor solution for the field equation derived from the more general 
Lagrangian 

m 

L = Co + Yl CkRn^R. (3.29) 

fe=0 
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Idea of Proof: Concerning fourth-order gravity this method was previously 
apphed e.g. by Barrow [64]. The de Sitter spacetime is an exact solution for 
the field equation, if and only if 2RG = (n + If we make the ansatz 

a{t) = 1 + (5{t) (3.30) 

we get the linearized field equation = /3 + n(3 for the Lagrangian . 
For the Lagrangian L — RD'^R we get the linearized field equation O'^R — 
(□'^i?) Q. For the characteristic polynomial we get a recursive formula such 
that the next order is received from the previous one by multiplying with 
■X ■ {x + n). We get the roots xi = —n — 1, X2 = —n (/c-fold), ^3 = 
(/c-fold) and x^^ = +1. Because of the last root the de Sitter spacetime is not 
an attractor solution. For the Lagrangian (3.29) we get the characteristic 
polynomial 



P{x) = x{x + n) 



Cq + Y, Ckx''~^{x + n)'''^{x - l)(a; + n + 1) 

k=l 



(3.31) 



for the linearized field equation. The transformation z = x"^ + nx + ^ 
gives a polynomial Q{z) which can be solved explicitly. So one can find 
those combinations of the coefficients such that the de Sitter spacetime is an 
attractor solution. 

It turned out that all the variants of the definition of an asymptotic de 
Sitter solution given in subsection 3.1 lead to the same class of solutions, i.e., 
the validity of the theorems written below does not depend on which of the 
variants of definition of an asymptotic dc Sitter spacetime listed is applied. 
For the 6th-order case we can summarize as follows: Let 

L ^ R'^ + CiRDR (3.32) 

and 

Le = R - - L (3.33) 

with length I > 0. Then the following statements are equivalent. 

1. The Newtonian limit of Le is well-behaved, and the potential consists 
of terms e""''^ /r with a > only. 
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2. The de Sitter spacetime with H = l/l is an attractor solution for L 
in the set of spatially flat Friedmann models, and this can already be seen 
from the linearized field equation. 

3. The coefficient ci > 0, and the graceful exit problem is solved for the 
quasi de Sitter phase H < l/l oi Le- 

4. P — 1"^ + l\ such that P Ci = l\ has a solution with Q <Iq <l\. 

5. < ci < P/A. 

A formulation which includes also the marginally well-behaved cases reads 
as follows: Let L and Le as in the previous result, then the following state- 
ments are equivalent: 

1. The Newtonian limit of Le is well-behaved, for the potential we 
allow 1/r and terms hke 

Pir) 

-^e"^' with q;>0 (3.34) 
r 

and a polynomial P. 

2. The de Sitter spacetime with H — l/l cannot be ruled out to be an 
attractor solution for L in the set of spatially fiat Friedmann models if one 
considers the linearized field equation only. 

3. Le is tachyonic-free. 

4. = /q + ll with P Ci = /q ll has a solution with < Iq < h- 

5. For the coefficients we have < ci < j. 

3.4 Discussion of no hair theorems 

In subsection 3.3 we have shown: The results of the Starobinsky model are 

structurally stable with respect to the addition of a sixth-order term ~ 

if the coefficients fufil certain inequalities. For the eighth-order case we got: 

For 

L ^ R'^ + CiRnR + C2ROOR, C2j^0 (3.35) 

and the case n > 1 the de Sitter spacetime with if = 1 is an attractor solution 
in the set of spatially fiat Friedmann models if and only if the following 
inequalities are fulfilled: 

< ci < < C2 < , (3.36) 

n + 1 [n + iy 
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and 

ci > -(n^ + n + 1)C2 + \l{n^ + 4n3 + An^^l + 4c2 (3.37) 

These inequalities define an open region in the Ci — C2— plane whose boundary 
contains the origin; and for the other boundary points the linearized equation 
does not suffice to decide the attractor property. 

This situation shall be called "semi-attractor" for simplicity. In a gen- 
eral context this notion is used to describe a situation where all Lyapunov 
coefficients have non-positive real parts and at least one of them is purely 
imaginary. 

In contrary to the 6th-order case, here we do not have a one-to-one 
correspondence, but a non-void open intersection with that parameter set 
having the Newtonian limit for Le wcll-bchavcd. 

To find out, whether another dc Sitter spacctime with an arbitrary Hubble 
parameter if > is an attractor solution for the eighth-order field equation 
following from the above Lagrangian, one should remember that H has the 
physical dimension of an inverted time, ci is a time squared, C2 is a time to 
power 4. So, we have to replace ci by ciH"^ and C2 by C2ii^ in the above 
dimensionless inequalities to get the correct conditions. Example: 

< CiH^ < ^ . 

Let us summarize: Here for a theory of gravity of order higher than four 
the Newtonian limit and the attractor property of the de Sitter spacetime are 
systematically compared. It should be noted that the details of the theory 
sensibly depend on the numerical values of the corresponding coefficients. So, 
no general overall result about this class of theories is ever to be expected. 

We have found out that for the class of theories considered here, one of 
the typical indicators of instability - cosmological runaway-solutions - need 
not to exist, even for an arbitrarily high order of the field equation. It is 
an additional satisfactory result that this takes place in the same range of 
parameters where the Newtonian limit is well behaved. 
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3.5 The Newtonian limit of 4th-order gravity 

Let us consider the gravitational theory defined by the Lagrangian 

Lg = {St^G)-^{r/2 + {aRijR'^ + PR^)f) . (3.38) 

G is Newton's constant, I a couphng length and a and (5 numerical param- 
eters. Rij and R are the Ricci tensor and its trace. Introducing the matter 
Lagrangian Lm and varying Lg + Lm one obtains the field equation 

Eij + aHij + pGij = SttG T^,- . (3.39) 

For q; = /3 = this reduces to General Relativity Theory. The explicit 
expressions Hij and Gij can be found in Stelle [59]. 

In a well-defined sense, the weak-field slow-motion limit of Einstein's 
theory is just Newton's theory. In the following, we consider the analogous 
problem for fourth-order gravity eqs. (3.38), (3.39). 

The slow- motion limit can be equivalently described as the limit c — > oo, 
where c is the velocity of light. In this sense we have to take the hmit G — > 
while G ■ c and I remain constants. Then the energy-momentum tensor T^j 
reduces to the rest mass density p: 

T,, = 6y^p, (3.40) 

x*^ — t being the time coordinate. The metric can be written as 

ds'^ = (1 - 2(j))df - (1 + 2ip) {dx^ + dy^ + dz'^) . (3.41) 

Now eqs. (3.40) and (3.41) will be inserted into eq. (3.39). In our approach, 
products and time derivatives of (p and ip can be neglected, i.e., 

i? = 4A^-2A0, where ^ = f,.. + f,yy + f,,z ■ (3.42) 

Further Rqq = -A0, Hqq = -2Ai?oo - Ai? and Goo = -4Ai?, where Z = 1. 

Then it holds: The vahdity of the 00-component and of the trace of eq. 
(3.39), 

i?oo - i?/2 + aHQo + /?Goo = SvrGp (3.43) 
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and 

-R - 4(a + 3/5) Ai? = SttCp , (3.44) 

imply the validity of the full eq. (3.39). 

Now, let us discuss eqs. (3.43) and (3.44) in more details: Eq. (3.43) 

reads 

-A0 - R/2 + a(2AA0 - Ai?) - 4/?Ai? = SttGp . (3.45) 
Subtracting one half of eq. (3.44) yields 

- A0 + 2q;AA0 + (a + 2/3) Ai? = AnGp . (3.46) 

For a + 2/3 = one obtains 

-(1 - 2q;A)A0 = AnGp (3.47) 

and then = is a solution of eqs. (3.43), (3.44). For all other cases the 
equations for and do not decouple immediately, but, to get equations 
comparable with Poisson's equation we apply A to eq. (3.44) and continue 
as follows. 

For q; + 3/3 = one gets from eq. (3.46) 

-(l-2aA)A(/) = 47rG(l + 2aA/3)p. (3.48) 

The A-operator applied to the source term in eq. (3.48) is only due to the 
application of A to the trace, the original equations (3.43), (3.44) contain 
only p itself. 

For a = one obtains similarly the equation 

-(l + 12/3A)A0 = 47rG'(l + 16/3A)p. (3.49) 

For all other cases - just the cases not yet covered by the literature - the 
ehmination of from the system (3.43), (3.44) gives rise to a sixth-order 
equation 

-(1 + 4(a + 3/3)A)(l - 2aA)A0 = At:G{1 + 2(3a + 8/3)A)p. (3.50) 

Fourth order gravity is motivated by quantum-gravity considerations and 
therefore, its long-range behaviour should be the same as in Newton's theory. 



34 



Therefore, the signs of the parameters a, (3 should be chosen to guarantee an 
exponentially vanishing and not an oscillating behaviour of the fourth-order 
terms: 

q;>0, q; + 3/3<0. (3.51) 

On the other hand, comparing parts of eq. (3.50) with the Proca equation it 
makes sense to define the masses 

m2 = (2q;)"^/^ and mo = (-4(q; + 2,(5))'^''^ . (3.52) 

Then eq. (3.51) prevents the masses of the spin 2 and spin gravitons to 
become imaginary. 

Now, inserting a delta source p — mS into eq. (3.50) one obtains for 
the same result as Stelle [59], 

(f) — mr~^{l + exp(— mor)/3 — 4exp(— m2r)/3) . (3.53) 

To obtain the metric completely one has also to calculate ■0- It reads 

ip — mr~^{l — exp(— mor)/3 — 2 cxp(— m2r)/3) . (3.54) 

For finite values mo and m2 these are both bounded functions, also for r — 0. 
In the limits a — > 0, i.e. m2 oo, and a + 3(3 ^ 0, i.e. mo — ^ oo. the terms 
with mo and m2 in eqs. (3.53) and (3.54) simply vanish. For these cases 
and ijj become unbounded as r — > 0. 

Inserting eqs. (3.53) and (3.54) into the metric (3.41), the behaviour of 
the geodesies shall be studied. First, for an estimation of the sign of the 
gravitational force we take a test particle at rest and look whether it starts 
falling towards the centre or not. The result is: for mo < 2m2, gravitation 
is always attractive, and for mo > 2m2 it is attractive for large but repelling 
for small distances. The intermediate case mo = 2m2, i.e., 3a + 8^ — 0, is 
already known to be a special one from eq. (3.50). 

Next, let us study the perihelion advance for distorted circle-like orbits. 
Besides the general relativistic perihelion advance, which vanishes in the 
Newtonian limit, we have an additional one of the following behaviour: for 
r — > and r — > oo it vanishes and for r ~ 1/mo and r ~ l/m2 it has local 
maxima, i.e., resonances. 
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Finally, it should be mentioned that the gravitational field of an extended 
body can be obtained by integrating eqs. (3.53) and (3.54). For a spherically 
symmetric body the far field is also of the type 



and the factors a and b carry information about the mass distribution inside 
the body. 

3.6 Higher-order gravity 

One-loop quantum corrections to the Einstein equation can be described by 
curvature-squared terms and lead to fourth order gravitational field equa- 
tions; their Newtonian limit is described by a potential "Newton -|- one 
Yukawa term". A Yukawa potential has the form exp(— r//)/r and was orig- 
inally used by Yukawa to describe the meson field. 

Higher-loop quantum corrections to the Einstein equation are expected 
to contain terms of the type RU^R in the Lagrangian, which leads to a 
gravitational field equation of order 2k + 4. Some preliminary results to 
this type of equations are already due to Buchdahl [38]. For A: = 1, the 
cosmological consequences of the corresponding sixth-order field equations 
are discussed in [66] and [75]. 

In the present chapter we deduce the Newtonian limit following from 
this higher order field equation. The Newtonian hmit of General Relativity 
Theory is the usual Newtonian theory. Prom the general structure of the 
linearized higher-order field equation one can expect that for this higher- 
order gravity the far field of the point mass in the Newtonian limit is the 
Newtonian potential plus a sum of different Yukawa terms. And just this 
form is that one discussed in connection with the fifth force. Here we are 
interested in the details of this connection between higher-order gravity and 
the lengths and coefficients in the corresponding Yukawa terms. 

Let us start with the Lagrangian 



mr 



-1 



1 -|- aexp(— mor) -|- 6exp(— 



(3.55) 





(3.56) 
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In our considerations we will assume that for the gravitational constant G and 
for the speed of light c it holds G = c = 1. This only means a special choice 
of units. In eq. (3.56), R denotes the curvature scalar, □ the D'Alembertian, 
and g the determinant of the metric. Consequently, the coefficient ak has the 
dimension "length to the power 2k + 2" . 

The starting point for the deduction of the field equation is the principle 
of minimal action. A necessary condition for it is the stationarity of the 
action: 

= 87rr^V^, (3.57) 



where T'^^ denotes the energy-momentum tensor. The explicit equations for 

P'' V^^-^ (3.58) 

arc given in [82]. Here we only need the linearized field equation. It reads, 
cf. [75] 

R ^ 

pij ^ j^ij g'^ + 2Y^ akig'^o''^^ R - O'^R''^] = SnT'^ , (3.59) 

2 fe=0 

and for the trace it holds: 

n-1 ^ 

9ij ■ P"' = + 2n ^ auig'^n^^^ R\ = 87rr . (3.60) 

^ fe=0 

n is the number of spatial dimensions; the most important application is of 
cause n — ?). Prom now on we put n = 3. 



3.7 The Newtonian limit in higher-order gravity 

The Newtonian limit is the weak-field static approximation. So we use the 
linearized field equation and insert a static metric and an energy-momentum 
tensor 

p>0 (3.61) 

into eq. (3.59). 



37 



Without proof we mention that the metric can be brought into spatially 
conformally fiat form, and so we may use 



9ij — Vij ~l~ fij ) 
Tjij = diag(l, -1, -1, -1) and 
fij = diag(-2$, -2*, -2*, -2*) . (3.62) 

Then the metric equals 

ds^ = (1 - 2^)dt'^ - (1 + 2^){dx^ + dy^ + dz^) , (3.63) 

where $ and \[' depend on x, y and z. Linearization means that the metric 
Qij has only a small difference to rjij; quadratic expressions in and its 
partial derivatives arc neglected. We especially consider the case of a point 
mass. In this case it holds: $ = $(r), ^ — ^('"), with 

r^^x^ + y^ + z^, (3.64) 

because of spherical symmetry and p = m6. Using these properties, we 
deduce the field equation and discuss the existence of solutions of the above 
mentioned type. 

At first we make some helpful general considerations: The functions $ 
and ^ are determined by eq. (3.59) for i = j = and the trace of eq. 
(3.59). If these two equations hold, then all other component-equations are 
automatically satisfied. For the 00-equation we need it!oo: 

Roo = -A$ . (3.65) 

Here, the Laplacian is given is as usual by 

g2 g2 g2 

dx"^ dy"^ dz'^ 

For the inverse metric we get 

g'^ = diag (1/(1 - 2$), -1/(1 + 2^r), -1/(1 + 2*), -1/(1 + 2*)) (3.67) 



38 



and 1/(1 — 2$) = 1 + 2$ + where is quadratic in $ and vanishes 

after hnearization. So we get 

gij ^ yfj _ fij . (3.68) 

In our coordinate system, /^-^ equals fij for all For the curvature scalar 
we get 

R = 2(2A* - A$) . (3.69) 

Moreover, we need expressions of the type o'^ R. OR is defined by OR = 
R-ijg^\ where ";" denotes the covariant derivative. Remarks: Because of 
linearization we may replace the covariant derivative with the partial one. 
So we get 

□'^ = ( 1)^= 2(-A*^+^$ + 2A*^+^*) (3.70) 
and after some calculus 

-Snp - A$ + A* . (3.71) 

We use eq. (3.71) to eliminate ^' from the system. So we get an equation 
relating $ and p — mS. 

-Att Ip + 8J2 aik(-l)*^A*^+^ p j = A$ + 6 ajfc(-l)*^A*^+2$ . (3.72) 

V k=0 / k=Q 

In spherical coordinates it holds 

A$ = ^$,, + (3.73) 

because $ depends on the radial coordinate r only. 

We apply the following lemma: In the sense of distributions it holds 

A ( -e-"/') = ^e-'-/' - AttS . (3.74) 

Vr / ri'^ 

Now we are ready to solve the whole problem. We assume 

TTL I ^ \ 

* = - 1 + E exp(-r//0 , k>0. (3.75) 



r 



i=0 
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Without loss of generality we may assume li ^ Ij for i ^ j. Then eq. (3.72) 
together with that lemma yield 



1 / n P 



k=0 i=0 \''» k=0 



-47r£ c, + 6£afc(-l) 



t 



k+l 



1=0 \ fe=0 

+2in±±±c,a,{-iy^'-^,A''^'S (3.76) 

k=0 j=k i=0 ''i 

where ti — ; therefore also ti ^ tj for i ^ j. This equation is equivalent 
to the system 



Eq = 1/3, 



i=0 



1=0 

0, s = 



■P 



(3.77) 
(3.78) 
(3.79) 



fc=0 



From eq. (3.79) we see that the values ti represent q + 1 different solutions 
of one polynomial. This polynominal has the degree p + 1. Therefore q < p. 
Now we use eqs. (3.77) and (3.78) . They can be written in matrix form 

as 

/ \ f Co \ / 1/3 \ 

l/to . . . l/tg 



Cl 







(3.80) 



V J 



Here, the first q + 1 rows form a regular matrix, the Vandermonde matrix. 
Therefore, we get 



l/^^ = EV/^f j = g + l,...p 

k=Q 



(3.81) 



with certain coefficients Xjk i.e., the remaining rows depend on the first q + 1 
ones. If Ajo 7^ then the system has no solution. So Xjo = for all g + 1 < 
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j < p. But ior q < p we would get 



(3.82) 



fe=i 



and this is a contradiction to the above stated regularity. Therefore p equals 
q. The polynomial in (14) may be written as 



/ 1 l/to-.-l/tg \ 
^ll/t,...lltl J 



( 



ao 



\ 



\ (-1)% / 



(3.83) 



This matrix is again a Vandermonde one, i.e., there exists always a unique 
solution (ao, . . . Op), which are the coefficients of the quantum corrections to 
the Einstein equation, such that the Newtonian limit of the corresponding 
gravitational field equation is a sum of Newtonian and Yukawa potential with 
prescribed lengths {{. A more explicit form of the solution is given in section 
3.9. 



3.8 Discussion of the weak- field limit 



Let us give some special examples of the deduced formulas of the Newtonian 
limit of the theory described by the Lagrangian (3.56). If all the vanish 
we get of course the usual Newton theory 



m 



$ = — , A$ = 
r 



-47r5. 



(3.84) 



$ and refer to the metric according to eq. (3.63). For p = we get for 
ao < 

m r 1 , ; = 1 

(3.85) 



m 

$ = — 
r 



1 I p-f-ZV-eao 

3 



cf. [59] and 



m 

* = — 

r 



1 



(3.86) 



For flo > no Newtonian limit exists. 

For p = 1, i.e., the theory following from sixth-order gravity 



(3.87) 
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we get, see [81] 



and 



where 



and 



m 
r 

m 



Co,i 



1 - CqC 



1 



r/lo 



-r/lo 



Cie 



cie 



-r/h 



-r/h 



ao 



6 2 J9a§ + 6ai 



Zo,i = Y^-3ao ± -^/Oag + 6ai . 



(3.88) 
(3.89) 
(3.90) 

(3.91) 



This result is similar in structure but has different coefficients as in fourth- 
order gravity with included square of the Weyl tensor in the Lagrangian. 

The Newtonian limit for the degenerated case lo — li can be obtained by 
a limiting procedure as follows: As we already know ao < 0, we can choose 
the length unit such that = —1/3. The limiting case da^ + 6ai — >■ may 
be expressed by ai = — 1/6 + e"^. After linearization in e we get: 

li^l±y/3/2eci^l/6±l/(6V6s) (3.92) 
and applying the limit £ ^ to the corresponding fields $ and ^ we get 



$ = m/r{l + (1/3 + r/6)e-'^} 
* = m/r{l - (1/3 + r/6)e-'^} . 



(3.93) 



For the general case p > 1, the potential is a complicated expression, but 
some properties are explicitly known, these hold also for p = 0, 1. One gets 



and 



$ = m/r I 1 + Q exp{—r/li) 

V 1=0 / 



'^ = m/r I 1 — Q exp(— r//j) 

V i=Q J 



(3.94) 



(3.95) 



where = 1/3; Yl, means Z^f=o ^« ^^<i ^^6, up to permutation of 
indices, uniquely determined by the Lagrangian. 
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There exist some inequalities between the coefficients a^, which must be 
fulfilled in order to get a physically acceptable Newtonian limit. By this 
phrase we mean that besides the above conditions, additionally the fields $ 
and ^ vanish for r — > oo and that the derivatives d^/dr and d^/dr behave 
like 0(l/r^). These inequalities express essentially the fact that the k are 
real, positive, and different from each other. The last of these three conditions 
can be weakened by allowing the q to be polynomials in r instead of being 
constants, cf. the example with p = 1 calculated above. 

The equality = 1/3 means that the gravitational potential is un- 
bounded and behaves, up to a factor 4/3, like the Newtonian potential for 
r ~ 0. The equation $ + \1' = 2m /r enables us to rewrite the metric as 



ds^ = (1 - 29) (1 - 2m/r)dt^ - (1 + 2m/r){dx^ + dy^ + dz^) 



(3.96) 



which is the conformally transformed linearized Schwarzschild metric with 
the conformal factor 1 — 2^, where 



e 



m 



E 



-r/k 



(3.97) 



can be expressed as functional of the curvature scalar, this is the hnearized 
version of the conformal transformation theorem. For an arbitrary matter 
configuration the gravitational potential can be obtained by the usual inte- 
gration procedure. 



3.9 A homogeneous sphere 

For general p and characteristic lengths /j fulfilling < /q < < 
write the Lagrangian as 



6 



(/^ + . . . + ll)R + {lili + lili . . . + li^,li)^R-r 



the coefficients in front of n^R in this formula read 



< Ip we 



(3.98) 



o<jo<h<-<ii<p "1=0 



(3.99) 
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Using this form of the Lagrangian, the gravitational potential of a point mass 
reads 



m 

$ = — 
r 



m 

•qj = — 
r 



(3.100) 
(3.101) 



For a homogeneous sphere of radius tq and mass m we get 



m 



1 



1 + - Ee"'^'' Ci(^o cosh(ro//i) - sinh(ro//i)) , (3.102) 

^0 i=0 



where 



(-ir^^nil-ir- 



(3.103) 
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